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Abstract 



In this paper we develop the Lagrangian and multisymplectic structures of the Heisenberg magnet (HM) model 
which are then used as the basis for geometric discretizations of HM. Despite a topological obstruction to the 
existence of a global Lagrangian density, a local variational formulation allows one to derive local conservation 
laws using a version of Nother's theorem from the formal variational calculus of Gelfand-Dikii. Using the 
local Lagrangian form we extend the method of Marsden, Patrick and Schkoller to derive local multisymplectic 
discretizations directly from the variational principle. We employ a version of the finite element method to 
discretize the space of sections of the trivial magnetic spin bundle N — M x S 2 over an appropriate space-time 
M. Since sections do not form a vector space, the usual FEM bases can be used only locally with coordinate 
£SJ i transformations intervening on element boundaries, and conservation properties are guaranteed only within an 
element. We discuss possible ways of circumventing this problem, including the use of a local version of the 
method of characteristics, non-polynomial FEM bases and Lie-group discretization methods. 
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1 Introduction 

The treatment of PDEs in invariant form finds its most natural language in the setting of jet bundles over an 
appropriate space-time. In this setting, governing equations of "motion" are prescribed by invariantly-defmed 
differential operators, and although calculations are ultimately done in terms of partial derivatives, the invariant 
meaning of all the operations is clear. Moreover, the invariant formulation and structure of the underlying fiber 
bundle frequently translate to a nontrivial geometry of the space of solutions, conservation properties and other 
"geometric" features. Although these geometric properties can be discovered in the coordinate formulation, 
conceptual clarity of invariant geometric derivations is frequently a great asset in itself. 

All of this is particularly so for Langrangian PDEs whose equations of motion are derived as extermality 
conditions on the action functional - an apparently coordinate-free form. Symplectic and multisymplectic struc- 
tures immediately follow by well-known procedures, as do conserved quantities corresponding to continuous 
symmetries (Nother's theorem). The situation is less fortunate for Hamiltonian systems as here a number of 
"symmetries" have been broken relative to the Lagrangian setting. A splitting of space-time relies upon the 
choosing a preferred time direction and a nonunique complementary space - a literal symmetry breaking as the 
full group of space-time transformations no longer preserves the form of the equations. Furthermore, the system 
is no longer specifed by prescribing a single quantity - apart from the energy functional a Poisson structure has 
to be specified independently 1 . 

Although it may appear from familiar examples of classical mechanics that the Lagrangian and Hamiltonian 
pictures are equivalent, this is not always so. In fact there is a local nondegeneracy condition on the Poisson 
bracket that must be satisfied before a local Lagrangian density can be reconstructed. In general, an additional 
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x In fact, there is a certain nonuniqueness associated with the choice of the canonical 1-form in Lagrangian theories living on jet 
bundles TV^' with r > 1 IH0, although this is not as bad as the number of choices that need to be made in the Hamiltonian setting. 



vanishing condition on a topological obstacle (i.e. exactness of the symplectic form) has to be satisfied before 
the variational picture is restored globally. Even so, a local Lagrangian formulation carries enough information 
to reconstruct the local multisymplectic geometry and Nother-type local conservation laws. A central result 
is the derivation of the local multisymplectic form for the Heisenberg Magnet model from its local Lagrangian 
formulation. 

Even when Lagrangian and Hamiltonian theories are equivalent, their discretizations may not be so. While 
the discrete variational principle yields analogs of symplectic and multisymplectic geometries under minimal 
conditions on the discretization, the Hamiltonian case calls for a sophisticated theory of symplectic or Poisson 
integrators to preserve the corresponding geometries. Therefore, the search for an equivalent Lagrangian formu- 
lation naturally arises when considering conservative discretizations. In this paper we discuss potential ways of 
constructing geometric discretizations from the local variational picture. 

In Section we briefly recall the picture of geometric PDEs as conditions on sections of an appropriate fiber 
bundle and their jets. The formal calculus of variations of Gelfand-Dikii is a very convenient computational 
tool for Lagrangian systems in this setting, and a source of many explicit formulas, including local conserved 
densities and a formal analog of Nother's theorem. Section[3]contains the derivation of the local Lagrangian form 
of HM, its conservation laws, and topological obstructions to the existence of a global description. The relation 
to Novikov's multivalued version of Morse theory as well as Lyusternik-Schnirelmann theory are discussed in this 
section. In Section 0] we discuss variational discretizations of HM using a version of the finite-clement method. 
An essential difficulty here is that the sections are constrained to lie on the unit 2-sphere, which is not satisfied 
by the standard FEM bases, and the resulting appearence of rational terms in the discrete Lagrangian. We 
discuss ways of avoiding this problem, including the use of discretizations based on Lie group methods to satisfy 
the sphere constraint in the most natural manner. 



2 Setting 

The geometric theory of differential operators and variational calculus is most naturally set within the framework 
of smooth fiber bundles N — ► M over some n— dimensional base space-time M and the corresponding jet bundles 

j\j(r) — > 21 . Differential operators and action functionals of the Lagrangian formalism act on the space 
of S{N) sections of the bundle, each section a : M — > N, tt o a = idM identifiable with its image E = cr(M), 
a copy of M horizontally embedded in N. Fixing a local trivialization N=\ oc M x F expresses A*" as a product 
of M with the fc-dimensional typical fiber F. A choice of local coordinates, x = (xi), t = l,...,n on the base 
and u = (uj),j = 1, . . . , r on the the fiber, defines a local coordinate system x x u : N -> R" +fe , making N 
locally diffeomorphic to the canonical projection n = ]R" +fe — -> R", while inducing a similar diffeomorphism 
n (r) s jjrc+fc x E n r _^ jpn for the N ( r ) with tlle a( i a pted coordinates ((xt), (u^)) > M < r2 on N(r) - The usual 

bundles A k TN and A k N have as their spaces of sections the k- vector fields V k (N) = S(TN) and the differential 
fc-forms D, k N = S(A k N). In addition, there is a distinguished subspace V{N,M) C V(N) = V 1 (N) of fiber- 
preserving fields, which generate local fiber bundle automosphisms, and the subspace V(N/M) = S(T(N/M)) C 
V(N, M) of sections of the vertical subbundle, which fix fibers. When restricted to E the local coordinates u 
define functions of x by 

m s : R" M A E A U k . 

Similarly, restriction of the fiber coordinates to E( r ) defines functions that encode the infinitesimal 

behavior of the section: (gjr)" ue(i) = u^ r) {x), \a\ < r. 

To the variational principle we admit local functionals A that to each a E S(N) assign the integral J A(cr) 

of an n- form A(cr) € Q, n M, such that the integrand A x (a) £ A™M at any x G M is determined by the r-jet ax 
of the section at that point 3 . In other words A is a differential operator S(N) — * Ct n M of order < r and by 
definition can be factored through a map A* induced on sections by a smooth bundle map A : — > A n M 
over M. 

S(N<- r >) A * ) il n M 

j:] [i ■ 

S{N) — A^ ]R 

2 Definitions for multi-indices are as usual: a — (ai, . . . , a n ), \a\ = X^i Qi J {w)" = air" 1 ' ' ' af - "™ ' = d" 1 ■ ■ ■ d" n , ti has 1 in 
the i-th place and zeros everywhere else. 
3 For a globally determined r. 



By a certain abuse of notation denned in we shall call A the symbol of the operator A and write 

A(<7) = J A*((7«), 

where in order to make sense of the integral we restrict integration to be over nice (see jH]) compact domains 
S M CM. 

The symbol A is identified with a section of the vector subbundle A°< n N^ = TT^*K n M C A™N^ of 
horizontal forms vanishing on the fibers of the ver tical subbundle T(N^/M). Then the action of A can be 
defined simply by integration of its symbol over holonomic 4 horizontal submanifolds E^- 1 = cj^^Em) C N^: 



A s = A(a)= / A. (1) 

In terms of this geometric description, variational calculus investigates the response of A to infinitesimal varia- 
tions of the holonomic submanifold E^. 

A generic variation Zs of E is represented by a (local) 1-parameter bundle automorphism group , with the 
induced (local) base automorphism group cf> M — tt o <^y, by restriction = Z\s of the infinitesimal generator 
Z e V(N, M). In coordinates we have 



alt 



tf N * = x( x ) + v(x, u ) = Y, x ^7r + Y, y i (*> u )jr- ( 2 ) 

t=o ox% i ou i 




The action of Z on f2AT naturally splits into the horizontal X and vertical Y parts. The horizontal part acts 
on functions Aq = C°°(N) by differentiation along E, i.e. via pushing forward to E, while the vertical part 
acts by differentiating along the vertical curves cr t (x) via the natural action of 4b sections cr* = 0^ o a o <fi~^: 

l \ 3 

j 3 3 \ i 

Geometrically, Y is the projection of Z onto the vertical subspace at T a r x \(N/M~) parallel to the horizontal 
subspace T X E = T x o~ • T X M, and X is the projection onto T X T,. 

The action of Z on the higher jets a^(x) is easily computed by interchanging the order of differentiation, 
for instance, 

j.i J1 j,z j.i 

Thus, the first jet components of Y^' are obtained from components of Y by horizontal differentiation along 
T^E, using the second jet information usj t i>i(x) to determine the horizontal subspace TE^ 1 ). The general rule 
follows recursively. 

Jet bundles allow one to abstract from particular sections by encoding their infinitesimal information point- 
wise, similar to the way tangent bundles allow one to manipulate vectors without referring to local integral 
curves. In general, will act on algebras of function / 6 A r = C°°(N^) by the horizontal holonomic lift 
di = Hor(^-) and has the property that it reduces to the total derivative ^fr upon restriction to a holonomic 
section: di — o~^ r K-^r-. Since, as seen in the example above, di uses the higher jet data, it acts on A r as a map to 
A r +\: on coordinates we have diu " = it^ Q+e, ' ) . Therefore the operators di must be regarded as derivations of the 
algebra A — inj lim/l r = (J A r of functions on the infinite jet space N = proj limiV^. All derivatives Der(.A), 
formally viewed as vector fields V(N) on N are given by the infinite series W = '^2 i Wi^r- + Y^j a Wj.a — fsjj 

and derivations di lie within the subspace V(N,M) of formally fiber preserving fields, thus generating the 
submodule Hor(iV) over K. Likewise the vertical holonomic lifts Ver(y) of fields on Y 6 V(N/M) generate 
Ver(iV) G V(N/M) C V(N, M) as a module over A: 

d^-^ + zZ^r^y Ver(Y) = ^Y^°eV(N,M), Y 3 e A . 



4 In contrast to iV = N, for r > not every horizontal submanifold E r C 7V (r ' is a lift of some E, since E^ r ' must satisfy the 
holonnmv constraint ( -§-) a iS, . (x) = ?i' a ^V.TV 




It is easily checked that Ver(iV) commutes with all <9j. 

Geometrically, fields in Hor(iV) reduce to V(M) and fields in Ver(iV) reduce to V(N/M) upon restriction 
to holonomic submanifolds C c N. In particular, the above Z, when regarded in V(N,M) uniquely 
splits into a "horizontal" and "vertical" part, 

^E^pE^-E 

i j j \ i 

Each part uniquely lifts to X = Hor(Z) E Hor(iV) and Y = Ver(Z) £ Ver(iV), and together define a unique lift 
Z = Hor(Z) + Ver(Z) eHor(iV) ®Yer(N)_ 

The direct sum Hor(iV) © Ver(N) C V(N, M), generates a split of the space of 1-forms ft 1 AT = inj limft 1 ^ 
and of the usual "de Rham" differential d into the vertical and horizontal parts: V^N = ft 1,0 jV ff^N and 
d = S + d respectively, where ft 1,0 7V = S A consists of forms vanishing upon contraction with horizontal fields, 
and ft 0,1 N = dA vanishing upon contraction with vertical vector fields. The geometric meaning of the horizontal 
and vertical differentials is that d represents differentiations along the horizontal subspace at each level E^and 
reduces to the differential in fl l M upon restriction to a holonomic section, while 5 vanishes, measuring the 
"virtual variations" transversal to holonomic T,( r f 

Operating formally we can deduce the properties of the calculus defined on ftiV by the pair of differentials 
H3- In particular, S and d anti-commute with one another, and on the generators we have 

dxi — dxi, dUj = uj a+c *^ 'dxi, df = difdxi, 

i i 

5uf = du\ a) - £ uf + ^d Xl , S Xl =0,6/ = ^ ~%Suf, 
and the module ftiV decomposes with respect to 5 and d 

n m N= M k ' l N= Y n k '°NAn°' l N, 5 : n k ' l N ^ fl k+u , d : ft*- 1 -» n k ' l+1 

k-\-l=m k-\-l—m 

into exterior subalgebras of vertical fc-forms VL k ' Q N — A\5 U \\ X \ •••> & u j" it an d horizontal Worms 

Q°' l ~NA[dx il ,. . .,dx it ] (in particular, ft']V = for I > n), with each subalgebra Q k ' l N filtered Qq' 1 N C • • • C 

£l k > l N C • • • by the subalgebras of k, Z-forms fi^'iV = AjSuj" 1 ', . . . , SUj" k \ dxi t , . . . , dxi,] of differential degree 

dc S = E s Kl <_£• _ _ _ _ 

The fields V(N) act on ftiV by Lie derivatives £-jyft = W(Q) — (di(W) +i(W) d)ft, with the only nonvanishing 

contractions being i(di) dxi — 1 and i( ® a) ) Su^ = 1. The easily observed commutation relations simplify 
calculations of £ Vcr (lv) and ^Hor(N) : 

5i(Hor(F)) +i(Hor(N)) S = di(Ver{N)) +i(Ver(iV)) d=0, 
Z(n) = (5i(Ver(Z)) +i(Ver(Z)) S)Cl+ (di(Ror(Z)) + i(Hor(Z)) d) ft. 

This forms the basis of formal variational calculus, developed in detail in here we sketch the necessary 

results. 

The main result is the formula for the variational derivative ©, which computes the action of Z in @ on 
the symbol of the action functional A € ft°'"iV( r ) via the holonomic lift: 

Z(A) = Z(A) = i(Vev(Z)) SA + dz(Hor(Z)) A, 

since i(Ver(N)) ft°'' = dft fe ' n = for any k, I. This shows in what sense the vertical part of Y acts on A as on a 
function, and the horizontal part X acts on A as on a horizontal n-form (see eq. (4.37)). In coordinates we 
have 

A = Adxi A ■ ■ • A dx n , A £ A, SA = ^ g u ( a ) A t j iXl a ... a dx n = 

dA 

— pr S(d a u 3 ) A dxi A • • • A dx n . 



By a formal differential analog of integration by parts we can split 5 A as follows: 

oA = oui — A ax i A • • • A ax n — 
5uf 



A • • • A dxi A • • • A da; n I € fiJ'"AT 8 dn (1,7l_1) iV. (3) 



The geometric meaning of the variational derivative operator is that it provides a formal adjoint to the operator 
of vertical lift on vectors Y G V(N/M), modulo the horizontal differential d: 

i(Ver(F)) <5A = i(Y) SA - i(Y) dQ = i(Y) SA a + d (i(Y) , (4) 

where 5Aq € flg' n N is the variational derivative. Integrating Q over some E^, since Y vanishes on the 
boundary of c>E, the adjoint meaning becomes exact and the equations of motion defining the extrema of A are 
given locally by 

5 A 

~r = 0- 5 

Suj 

Define the multisymplectic form Q and the vertical part Xi of di by 

n = 5e e n 2 ' n -\ nW = <5e w e n 2 ' , x t = a,- - A = V u[ a+u) e v(n/m). (6) 

OXi J 

On any holonomic E( r ) satisfying (JjJJ we have the equivalent covariant Hamiltonian system, which, when A does 
not depend on the base coordinates x explicitly, has the form 5 

^i(Xi)n® =-8H, H = J2 i ( X i) @(j) ~ Ae A ( ? ) 

i 3 

obtained by a version of the Legendre transform. Among other things, it implies the preservation of the multi- 
symplectic form 

J2Xi(Q^) = -6 2 n = Q, (8) 

i 

as well as the local conservation laws (formal Nother's theorem): 

J2 d k T jk = 0, Tjj =H-J2 i&k) W , T jk = i{X 3 ) 9« j ? k. (9) 

k k^j 

Later we will have a chance to use the momenta defined as [3] 

= = V ^ H l^Tn (10) 

in particular, P Ji(0) = 

The connection with the usual time formalism is done as follows |3J. Fix a space-time split of M = M™ -1 x M 1 
with a compatible coordinate system (x\, . . . ,x n —i,t = x n ). For any section S fix the time slice E(i') over the 
hyperplane M n ~ x x t defined by x n = t. For the holonomic lift T,^ r '(t) of any such slice, we define the symplectic 
form lo, the corresponding primitive 1-form 8 and the Hamiltonian functional H with the corresponding densities 

&(o-)= / u, 0(a) = 9, H(tr)=/ T nn dx\ A ■ • • A dx„_i, (11) 

JSW(i) >/EW(i) JEW(() V v ' 

H 

provided the Lagrangian A is regular |SJ|S]. Then the Hamiltonian evolutionary equations are given as usual by 

<5H = i(X t )£, 

and the Lagrangian density is recovered by a version of the inverse Legendre transform: 

A = (-n + i(X t )6)Adx n . (12) 



5 Whinh dpnfinds on thp. chnirc of trivia.liza.tion. 



3 Local Lagrangian of HM System 



In this section we derive a local Lagrangian form of the action functional of the Heisenberg magnet (HM) model 
since the symplectic form corresponding to the Poisson bracket of HM is not exact and a global Lagrangian does 
not exist. Throughout this section we explicitly work in coordinates. 

The HM model is a Hamiltonian system defined on the product fiber bundle N — M x R 3 , where the space- 
time base is simply M — R 2 with coordinates (x%, x 2 ) on the base and S = (Si, S 2 , S3) on the fiber (e.g. [S])- 
A solution of the system is a section £ S N represented by the restriction of the fiber coordinates S : M — > R 3 
and satisfying the following equation 

9 2 S a — 2^ e abcSbdfS c = 0. 
b,c 

We now reconstruct the Hamiltonian structure of HM in two different coordinate systems on the fiber. 



3.1 Hamiltonian form of HM 

After the space-time split M = M n ~ x x M 1 = R x R with coordinates x = Xi and t = X2, we have the 
instantaneous state of the system E(i) = a(t)(M n ~ 1 ) at time t £ M 1 given by a section a(t) of the instantaneous 
bundle N(t) = M"" 1 x R 3 -> M ra_1 and represented by a function S(t) : M"" 1 R 3 . THe notation from 
Sectional is modified accordingly: A r (t),A(t), fUVW(t) etc. The Hamiltonian functional of the system is 



H(<r(t)) = / i {SxJ + S 2 . x 2 + S 3 J) dx, 



prescribed by the horizontal 1-form H = Hdx = \ (S hx 2 + S 2 , 2 + S 3 , x 2 ) dx e Q^N^it) e Q^Nit) with the 
density H S A\(t). There is a local Poisson bracket on each fiber defined in coordinates by 

{S a , S b } = * a , b (S) = ~Y, tabcSc, (13) 

c 

which can be viewed as a section over N(t) of the exterior power of the vertical bundle "J : N(t) — > /\ 2 T(N(t) /M). 
Since we are interested in a local description, we shall avoid the development of the usual infinite-dimensional 
formalism via the introduction of a Poisson bracket as a local bilinear operator ^(S(x))S(x — y) on the loop 
algebra of state variations (see |S]), and instead immediately write the equations of motion in the local form for 
any point x € M (t) : 

S t (x) = {H, S}(x) = ^a,h(S(x))^-(x)^-(x) = S(x) x S xx (x), (14) 

where as functions in A$(t) the coordinate functions S have the natural variational derivatives as defined in 
Section H f = E t M S i 



To restore independence of the fiber coordinates in (|14|) we can state that the vector product in (|14f> is 
defined with respect to the standard metric on R 3 . There is a more natural way, however: the bracket (| 1 3|) is 
the Lie-Poisson bracket on the space so(3)* = R 3 dual to the Lie algebra so(3) [5], and therefore variational 
derivatives of functions in A(t) = n°-°N{t) lie in Slo' 7V(t) = A(t) <x> T*so(3)* ^ A(t) ® so(3)*, on which the 
bivector field ^ £ A(t) ® A 2 Tso(3)* acts naturally and independently of any coordinate system. 

The bracket (|13(l is nondegenerate on the leaves of a 2— dimensional foliation of R 3 consisting of spheres about 
the origin, orbits of the coadjoint action of SO(3) 0. On each such leaf a symplectic form dual to the bracket, 
the Kirillov form, is defined. Since, as is easily verified, the dynamics of Q14[l preserves the leaves (\S\ 2 = 0), the 
system can be considered on the bundle N(t) — M n_1 x S 2 with the unit sphere as the fiber, provided the initial 
conditions are chosen as a section of this new bundle. It is well-known that the symplectic form dual to 'J is the 
area form on the sphere, which can be written down as soon as local coordinates have been chosen on S 2 . For 
instance, the unit sphere with a deleted north pole can be parametrized by the inverse stereographic projection 
from the pole (0,0, 1). Following [S] we introduce the following transformation to complex coordinates: 

w = wi + iw 2 , S = Si + iS 2 , Q — S 3 , R = 1 + ^ > 0, S = (S, S, Q) T , W = (w, w) T , 



the correspondence is then established as follows: 



2w _w _ 2w _ w _ \w\ 2 — 1_R— 1 
= 1 + H 2 = R' = 1 + H 2 = R' = M 2 + 1 = R ' 

H 2 = ^, R=t 1 t,, w = -^- = SR,w = SR. 



(15) 



The Poisson brackets of the modified coordinate functions S = (S, S, Q) and the Hamiltonian density 



{S, S}==2iQ, {S,Q} = -iS, {S,Q} = iS, H = -(\S X \ 2 + Q 



2 ) 



(16) 



generate the equations of motion equivalent to l|14|) . The coordinate transformations (|15|l generate the bracket 
of the W-coordinates: 

{w, w} = R 2 {S, S} + SR 3 {S, Q} + SR 3 {Q, S} = -2iR 2 . 
This bracket is clearly non-degenerate and defines a symplectic (Kirillov) form 

Sw A Sw 2Sw A Sw 



US = 



2iR 2 i(l + H 2 ) 2 ' 



(17) 



w — , , SzASz . . 

151 reduces to the canonical form : . Since exterior differentiation commutes 



which after scaling z 

with pullback, we obtain a primitive for us by pulling back a primitive 9 of the canonical form: 



zSz — zSz 1 
2i ~~ 2i 



w w\w\ 



R AR 2 



Sw 



The density in l|16f> transforms to 



4i? 2 



H 



■Sw — 



R 4R 2 



r w\wr r _ 
Sw + - —5 Sw 
AR 2 



wSw — wSw 
2iR ' 



w-w. 



2R 2 ' 



which generates equations of motion in the new bracket 



w t = 1- 



-Rw xx + ww x 
R 



-w xx - \w\ 2 w xx + 2ww x 
2R 



-iw xx + t- 



.ww x 



R 



(18) 



(19) 



3.2 Lagrangian and Multisymplectic formulation of HM 

We now have all the elements necessary to reconstruct the Lagrangian and multisymplectic form of the system 
in the complex plane. With respect to the symplectic structure (|17fl the equations of motion (|19f) are 

From l|12fl . letting 9 be defined as in i|3.1|) . X t = £>t = w t-§^ + ™fg=, and H = Hdx using l|18() . one obtains 



(H + i(X t ) 9) dx A dt 



W X W X WWt — ww t 



2R 2 



2iR 



dx A dt . 



(20) 



The multisymplectic form now follows from Q . Computing the variation of the Lagrangian one obtains 

Sw f\dx f\dt 



. ww x w x 

6A=< ^R^ + dx \2R 2 J 2iR 



( W x \ Wt WWWt — W 2 Wt „ / 

I 1 1 o t 



MR 2 



2iR 



+ 



'JW X W X / W x \ W t W 2 W t - WWWt 

2R 3 X \2R 2 ) 2iR AiR? 



(——) f Sw A dx A dt 
\2iRJ 



W x W x W W 

d S TTT^Sw Adt+ -—^Sw /\dt+ t~-^Sw A dx - -—=Sw /\dx>, (21) 
2R Z 2R l 2iR 2iR 



yielding the fundamental differential forms 
9 



w x r 1 w x f_ , W r , W r _ 

— -Sw A at -Sw A at Sw A dx H Sw A dx, 

2R 2 2R 2 2iR 2iR 



(22) 



_ WW X — WW X - r _ lc-r_, 

= S& — 3 SwASwAdt + —pSwASw x Adt 



2R2 



5u> A (5w x A c?i 



2?,7? 2 



(5wA<5wAdx, (23) 



and, by using Q, the covariant Hamiltonian 



n = 



w x w x ww t - ww t 



2R 2 



2iR 



dx A dt 



wwt — WWt w x w x 



2iR 



R 2 



dxhdt = ^dxAdt = -H. (24) 

2R 2 v ' 



Comparing the coefficients of Sw A dx A dt and Sw A dx A dt in 



S nj WW X W X . , 

<W = - „ n dw AdxAdt 



ww x w x 



2R 3 



„ -Sw Adx Adt — — %rSw x AdxAdt — — %r5w x AdxAdt, 
2R 3 2 R 2 2R 2 



to those in 



j A ./5 \r> + WW X W X - WW X 2 Rw xx + WW X W X - WW X 2 

dx A i(o x ) il = ^-^3 ow AdxAdt + ow AdxAdt 



2R 3 

-^^5w x Adx Adt — -^^Sw x Adx Adt, dtAi(d t )Sl 



2i? 3 



* -Sw A dx A dt — t„ Sw Adx A dt, 



2iR 2 2iR 2 

we obtain the equations of motion 119|) , while the coefficients of Sw x Adx Adt and Sw x Adx Adt yield the tauto- 
logical identities 

w x w x w x w x 

~2R 2 ^~2R 21 ^W^^^W 2 ' 
These can be used to define momenta as in <|10[) : 



SA dA 



Sw x dw x 



OH 

dw x 



w x _ _ SA _ dA 

~w?' p = Fw) ~m x ~m: 



dH w x , . 

aw x 2R Z 



so that in new fiber coordinates (w, w, p, p, Wt, Wt) we have 



H = -2R 2 ppdxAdt = -(1 + \w\ 2 )\p\ 2 dxAdt, Q = -SwASpAdt - Sw A Sp A dt + ——rSw A Sw A dx, (26) 

2iR 2 

and the components of 10 take the form 

dx A i(d x ) = w x Sp A dx A dt — p x 5w A dx A dt + w x 8p Adx Adt — p x Sw Adx Adt, 

dt A i(d t ) = „ W * Sw A dx Adt— W * Sw Adx A dt, 
2iR 2 2iR 2 

SH = -2Rw\p\ 2 Sw AdxAdt- 2Rw\p\ 2 Sw AdxAdt- 2R 2 p5p A dx Adt - 2R 2 pSp Adx A dt, 
yielding the equations 



-2R 2 p = w x , -2R 2 p = w x , -2Rw\p\ 2 = - 



2iR 2 



Px, 



-2Rw\p\ 2 



2iR 2 



(27) 



which are equivalent to l|19|l . This is the local multisymplectic formulation of HM. 

To relate our formalism to £Q we introduce two matrices K and M and a coordinate vector z = (w, w, p, p) 
(since Wt, Wt never show up explicitly, we can get away with fewer fiber coordinates): 



K 



[0 -1 0\ 

-1 

10 

\0 1 / 



M 



2iR 2 



/0 -1 0\ 

10 



\o o o oy 



(28) 



If Q = £l t dt + Q x dx, where Of, Q x S Q( 2 >°) are fiber 2-forms, then K and M are simply matrices for f2 t and — fi a 
in z-coordinates. Indeed, 

dx A i(d x ) = -i(d x ) Q, t dx Adt = -z x K • Sz A dx A dt, dt A i(d t ) = i(B t ) VL x dx Adt = -zj M • Sz A dx A dt, 



and the canonical equations (J2J) reduce to 



K • z x + M ■ z t = V Z H, 



where V z 7i is a column vector of partial derivatives. 

Finally, the conservation laws © yield the usual energy (H) and momentum (P) conservation laws: 



_ w x w x _ w t w x + w t w x 

J-tt — ~ r, — — til J-tx — —- 



T x t — 



2R 2 
ww x — ww x 
2~iR 



2R 2 

W X W X WWt — WWt 



2R2 



2iR 



(29) 

(30) 
(31) 



It is important that the definition of these quantities ultimately relies on the choice of the primitive form 9, 
which is defined only locally. In general, we have no global definition for the conserved densities Tjk nor for the 
Lagrangian density itself. The functionals defined by these densities are known as multivalued junctionals since 
changes in local choices of 9 result in differences by a multiple of a topological term, much like the multivalued 
function lnz. This phenomenon was investigated, among others, in [51 where [5] sets a general task of 
relating the properties of critical points of such functionals to the topology of the underlying infinite-dimensional 
manifold of sections, extending the classical Morse theory, and, similarly, in the finite-dimensional case F, 
replacing functions of the Morse theory by nontrivial elements H 1 (F, R), i.e. close but non-exact forms. In our 
situation the finite-dimensional case (9) fully determines the infinite-dimensional counterpart. The essential idea 
of |5| is to find a covering manifold on which the 1-form becomes exact, defining a global function, to which the 
classical Morse theory applies. At the moment it is not clear to us what covering of the sphere S 2 will trivialize 
9. Such a covering is unlikely to be obtained by reduction of the symmetry groups 50(3) and SU(2) since, being 
3-dimcnsional they do not carry a natural symplectic structure themselves, and being compact they do not have 
factor-spaces F with a nontrivial H 2 (F,WV). The necessary covering must be similar to an infinitely-sheeted 
surface of lnz, except covering the sphere completely (no branch points excluded). Further, the classical theory 
due to Lyusternic-Shnirelmann extends the Morse theory in the infinite-dimensional case to allow degenerate 
functionals, and relates the numbers of its critical points of different index to a topological invariant called the 
Lyusternik-Shnirelmann category (see the recent work by Walter Craig on water waves). This theory does require 
a globally-defined Lagrangian function, which the Novikov-Morse theory promises to provide. Given this, we 
should be able to estimate the number of solutions of a given kind (say, periodic ones) of HM on a given domain 
(e.g., with periodic boundary conditions) including the higher-dimensional spatial cases. Similar considerations 
apply to linearizations of HM near fixed points or periodic solutions, promising to be a tool of stability analysis. 
These investigations will be continued elsewhere. 



4 HM discretizations 

Given a PDE in variational form there are approaches to its discretizations that more or less automatically 
produce discrete analogs of the preservation of the multisymplectic form (JHJ, and frequently result in superior 
discrete conservation laws (Cf. © which suggests a discrete analog of Nother's theorem). While [H] developed 
a finite-difference approach to variational PDEs, we believe a more natural approach is to rely on the finite 
element method (FEM) involving space-time meshes. One of the reasons is that the equation already possesses 
a natural variational form, unlike many equations treated by a version of Galerkin's method. Here we outline 
our FEM-based approach to discretization of the HM model. 

4.1 Finite element method for HM 

Using H2()(l the continuous action functional is (with C = A as defined in QJ) 

r /"a f ( ■ wW t ~ Ww t w xW x \ 

c= J A= j [ i 2R ~wr* A ^ 

Introducing a rectangular space-time mesh M on M with elements {e}, we take a basis of piecewise linear on 
e functions P 1 (more generally they can be piecewise polynomial) 12] which, on the canonical square element 
with coordinates {xi,rj) 6 [0, 1] x [0, 1], have the form 

^ = \(l-t)(l-V), 02 = ^(1+0(1-'/), <t*=\(l+t)(l+V), 04 = ^(1-0(1+'/)- 

Then the approximations w = ^2 V w v 4> v are defined by collocation at element nodes (1 per element vertex, 4 per 
element), and for polynomial terms of C they reduce to polynomials in the nodal values w v , while the rational 
terms have no canonical representation. First, we consider the product approximation J2J j=r = r v (f> v = -^(f> v , and 



obtain the following discrete action functional: 

L = \j E ir u cj> u w v uf {cj> v dtfe ~ d t cf> v <fc) - £ (f"rV#^^^ d x fc) 



u.u.v.v 



dx A dt 



a, a, 6—1 



X! X! KK^^Ji^eMdt^e.a-dt^e^^e.ld^bdxhdt- ^ « 

v a,a,b,b— 1 



e' e 



i J d x (f) e . a d x (f) e .a(f> e , b <j) e j dxAdt 



H' 



a,a,b,b 



composed of element junctionals L e : 



4 4 

L = ^L e (t«e,le), L e {w e ,W e ) = i <<rlK,a,b- < <^ H l^ 

e a,a,6— 1 



a, a, 6. 6—1 



A e (u) e ,U) e ) 



called element Lagrangians. The coefficients A e a - b and H^_ b - are independent of the element e taken to be a 
rectangle with sides h x , h t , and which can be computed by reduction to the canonical integrals: 

A< a,a.b = Aa.a.b = \~ (j>a fyfe ~ d V 4>a <fe) 4>b ^ A df], ffj^j = = ~ ^ 0{<M& (j^ d£ A df], 

employing the obvious symmetries: 

Aj,a,fc = -^a,a,b,b G ]R ' ^0,0,6,5 = H a.a.bl = H a.a,b,b G IR - 

Being polynomials of degree < 8 in £, 77 these can be computed exactly (to round-off) using 5-point Gaussian 
quadrature, and as in any FEM code are never derived explicitely but assembled at run-time. While H is a full 
4-th order tensor, the 3-rd order tensor A has the following sparsity structure: 



f° 





X 


x\ 








X 


X 


X 


X 








\x 


X 





V 



As an alternative to the product approximation, we can discretize the rational terms by using their element 
averages. The element average of a function is defined to be the average of its approximation in the space 
spanned by {4> v }- 



4 4 

|C| Je v 11 Je a=l o=l 11 Je a 



1 = h x h t 



/° denotes the nodal value of / at the a-th vertex of clement e. Computing the element averages easily yields: 
<i>a = j,a = 1,.. .,4; then 



fe = l (f 1 e +f 2 e +f! + ft)- 



Defining element moment matrices K^>: 



(32) 



we have 



£(0) 





/4 


2 


2 


1\ 


h x h t 


2 


4 


f 
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36 


2 


1 


4 


2 




V 


2 


2 


V 



12 



(-2 -f 2 1\ 

-1-2 1 2 

-1-2 1 2 

\-2 -f 2 1/ 



6/l T 



/ 2 -2 1 -1\ 

-2 2 -f f 

1-1 2-2 

f -2 2) 



■ (33) 



Now substituting approximations w = ^2 V w v <ft v and R e = R(w e ) into C (R is defined in the element interior 
only, which is enough for the integral to make sense), we obtain another discrete action functional: 



E 



W V w"(j) v 4>v . W V W V (4> v (j) 1 j t t - 4>v,t4>v) 



2R 2 



2iR 



dxf\dt = ^2 E ( ~ 



«fA2)_ W >e ( r>W frtth 



2R 2 



2iR e 



L e (w,w) 



Noting that 

£(1) - (RW f = 

and letting K = we can write 

L = }^L e (w,w), L e (w,w) 



4 



I 
-I 



E 



2R P 



R, 



(34) 



The element Lagrangian L e corresponds to La in 8 and represents one element's contribution to the action. 
For a regular mesh the discrete Lagrangian is invariant under a mesh shift and as can be easily seen from 
l|34fl . then the element Lagrangian is defined in terms of the canonical element Lagrangian L: 



L e (w,w) = L(wl, wj, w^, wf), L{w l , w 2 , w 3 , w 4 ) = 



2R 



K„ 



R 



(35) 



where (w\, w 2 ,w 3 ,w 4 ) denote the four nodal points at a particular element, and R is defined using the average 
of the four values w 1 , w 2 , w 3 , w 4 . Defining J(w, w) and K(w, w) by their matrices 

J(w, W) — W ■ J ■ W = W a J a} aWa, K(w, W) = W • K ■ W = ^ W a K a ,aW~a, 



we obtain 

L(w, w) — -LrJ(w,w) \—K(w, w), 

y ' 2R y ' 2R 2 K 1 2R 



{w 1 w 3 — w 1 w 3 ) (w 2 w a — w 2 w 4 )~ [ 



1 th 



\2\w 



1 |2 



2\w 



2 1 2 



2\w 



3|2 



2\w 



4 I 2 



2R 2 Qh x 

- 2{w 1 w 2 + w'i 2 ) + (wV 3 + w 1 w 3 ) - (u/uJ 4 + W 1 w 4 ) ~ (w 2 w 3 + w 2 w 3 ) 

+ {w 2 w 4 + w 2 w 4 ) - 2{w 3 w 4 + w 3 w 4 )] . (36) 

The discrete Euler-Lagrange field equations (DELF, [H]) are obtained by differentiating the discrete action, 
which depends on a particular nodal value w v only through the element Lagrangians corresponding to four 
elements e", a = 1, ... ,4 containing v as one of their vertices: L e i , L e 2 , L e 3 , L e 4 . DELF equations are written 
simply as 

dL 



dw v 



-^=0 

' dw v ' dw a 8 



1^ 

w, 



determined using the canonical element derivative 



dL 



—~Ja ■ W 



dw a 2R 
4.2 Alternative approaches 



— — J(w, w) —K a ■ w H —K(w, w). 

WR 2 y ' 2R 2 8R 3 V ' 



Preliminary numerical experiments indicate that the FEM-based methods developed above can be unstable, 
apparently due to a relatively poor approximation of the rational terms in the polynomial bases. Here we sketch 
alternative approaches that may help avoid these difficulties. 

First, consider a local quasi-linear form of the HM equations of motion. For this construction we rewrite the 
Lagrangian for HM and the equations of motion in the real coordinates w = a + ib: 



W t W - WWt w x w x 



2iR 



2R 2 



dx A dt 



Im 



( ww t 



Re 



w x w x 
2R 2 



dx Adt = 
( 

ab t — a t b a 
R 



b 2 

X ' X 

~2R 2 ~ 



dx A dt. 



Applying the Euler-Lagrange operator explicitly we obtain 

SA b t fb\ abt-<Hb 2b t b , a 2 b t a t ab 2b t (a 2 + b 2 )b t b t 

77= d+ 9 ' d " a m = -5"- m (aat + bbt) 



5a R \R J R 2 R R 2 + R 2 R R 2 R 2 ' 

SH f a ^\, a l + b l a xx a x al + bl a xx a (b 2 x - a 2 x ) - 2a x bb x 

^ =d *{W +a -R^ = lp- 2 & {aa * + bb * ) + a -R^ = -R 2 - + W • 

Since A is anti-symmetric in a and b while H is symmetric in the same variables, we easily obtain the corre- 

spondmg expressions tor — — and -r— : 

06 oo 

— - h. &H _b xx b{a 2 x - b 2 x ) - 2aa x b x 
Sb ~ i? 2 ' Sb ~ R 2 i? 3 
Thus, the Euler-Lagrange equations are 

SA _ bt a xx a{a 2 x -b 2 x ) + 2a x bb x SA _ a t b xx b(b 2 x ~ a 2 x ) + 2aa x b x _ 

5a ~ R 2 R 2 i? 3 ' Sb ~ R 2 R 2 i? 3 ' [ ' 

which arc equivalent to l|19|l or 

dta + d x b x+ h ^- a \ +2aa * b * =0, dtb-8 xax+ < a *- b °l +2a * bb * =0. (38) 
R R 

This equation has the quasi-linear and even the semi-linear form and is amenable to the method of characteristics. 

Another approach, relying on the construction of a special class of section bases spaces relies on Lie-group 
methods. Limitations of space allow us only to sketch the approach. The group SO (3) acts transitively on S 2 , 
therefore on each element e with canonical local coordinates e have S (x, t) = g(x, t)- So, where g € SO (3) and the 
action should be coadjoint (unlike the standard linear representation, found in the numerical literature), since it 
interacts naturally with the commutator. The group near the identity must be parameterized by its Lie algebra 
= H 3 , so that g(x,t) — exp(xA(x, t) + tB(x,t)) where A and B are polynomial functions of x,t with values 
in so(3) and exp is an approximation to the exponential map (e.g. the Caley transform (/ + \A){I — ^A)^ 1 ). 
Using a convenient choice of A and B we should be able to reduce the spin length constraint (orbit throught 
So diffcomorphic to a the coset space by a close subgroup of rotations around So), but is linearization (ideal in 
so(3)). The spaces of sections (no longer linear) generated by A and B in a certain smoothness class replace the 
traditional FEM bases. The details are the subject of a forthcoming publication. 

In this paper we have established the local Lagrangian and multisymplectic structure of the Heisenberg 
magnet model, and have shown how the powerful formalism of the variational calculus can naturally lead to the 
conservation properties and suggest natural discretizations. 
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